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Abstract. We consider entanglement in the ground state of the XY model on 
infinite chain. We use von Neumann entropy of a sub-system as a measure of 
entanglement. The entropy of a large block of neighboring spins approaches a constant 
as the size of the block increases. We evaluate this limiting entropy as a function of 
anisotropy and transverse magnetic field. We use integrable Fredholm operators and 
the Riemann-Hilbert approach. The entropy reaches minimum at highly ordered 
states but increases boundlessly at phase transitions. 
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1. Introduction 



Entanglement is a primary resource for quantum computation and information 
processing [221 I2ni HI HI IHl IH]- It is necessary for quantum control. It shows how 
much quantum effects we can use to control one system by another. Stable and large 
scale entanglement is necessary for scalability of quantum computation [23 123 121] • 
Entropy of a subsystem as a measure of entanglement was discovered in X]- Essential 
progress is achieved in understanding of entanglement in various quantum systems [H 

la la 1101 iini ini ini iiHi iia 123 123 123 iisi 111 1201 

XY model in a transverse magnetic field was studied from the point of view of 
quantum information in [HI HI IS1 inH] ■ In this paper we evaluate the entropy of a block 
of L neighboring spins in the ground state of the XY model in the limit L — oo 
analytically. Our approach uses representation of 1^ and the Riemann-Hilbert method 
of the theory of integrable Fredholm operators. The final answer is given in terms of 
the elliptic functions and is presented in Eq. (pUj) and ()29|) below. The Hamiltonian of 
the XY model can be written as 

oo 

n = - E (1 + 7)«+i + (1 - iX^l+i + ha:, (1) 

n=— oo 

Here 7 is the anisotropy parameter; a^, and are the Pauli matrices and h is the 
magnetic field. The model was solved in ^21 Hll HHl (HH] • The methods of Toeplitz 
determinants, as well as the techniques based on integrable Fredholm operators, were 
used for the evaluation of some correlation functions, see [i^ t IHI H HH j HH j liH]. 

We consider the ground state \GS) of the model. We evaluate the entropy of 
a sub-system of the ground state. We shall calculate the entropy of a block of L 
neighboring spins, it measures the entanglement between the block and the rest of the 
chain We treat the whole chain as a binary system \GS) = \A^B). We denote this 
block of L neighboring spins by subsystem A and the rest of the chain by subsystem 
B. The density matrix of the ground state can be denoted by pab = \GS){GS\. The 
density matrix of subsystem A is p^i = TrBipAs)- The von Neumann entropy S{pa) 
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of subsystem A can be represented as follows: 



S{pa) = -Tva^Pa'^^Pa) 



[XI 



(2) 



This entropy also defines the dimension of the Hilbert space of states of subsystem A. 

A set of Majorana operators were used in with self-correlations described by 
the following matrix: 



B, 



Ho n_i 

Hi Ho 

n^-i . . . 



Hi 



no 



Here 



ZTC Jo 



9ie) 



and g{9) 



cos 6 — i'y sin 6 — h/2 



(3) 



cos 6' — i7sin^ — h/2\ 
One can use an orthogonal matrix V to transform B/, to a canonical form: 

l\ 

-1 

The real numbers — 1 < z/,„ < 1 play an important role. We shall call them eigenvalues 
The entropy of a block of L neighboring spins was represented in ^ as 



(4) 



S{pa) = E H{y. 

m=l 



(5) 



with 



H{v) 



1 + z/, 1 + z/ 1-z/, 1-z/ 
In In ■ 



(6) 



2 2 2 2 

In order to calculate the asymptotic form of the entropy it is not convenient to 
use formula (jH) and © directly. Following the idea we have already used in jlOj, let 
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us introduce: 



Block Toeplitz Determinants 



Bl{X) = iXh - Bl, Dl{\) = det Bi(A) 



(7) 



and 



e(x, u) 



X + V X + V X — u X — V 
m m ■ 



2 2 2 2 

Here II is the identity matrix of the size 2L. By definition, we have H{y) = e(l, u) 
and 

L 



^L(A) = (-i)^n(^ 



2 2 \ 
V 

ml 



(9) 



m=l 



With the help of the Cauchy residue theorem we rewrite formula (0) in the following 
form: 

S{pa)= lim -^i dAe(l + e,A)-^lnDi(A) * (10) 
e^o+ 47ri Jv dA 

Here the contour V is depicted in Fig [2 it encircles all zeros of Dl{X). We also notice 




Figure 1. Contours T' (smaller one) and T (larger one). Bold lines (— cx),— 1 — e) 
and (1 + £,00) are the cuts of ttie integrand e(l + e, A). Zeros of Dl(X) (Eq. 0) are 
located on ttie bold line (—1,1). The arrows indicate the directions of integrations, 
and r and R are the radius of the circles. % 
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that B/^(A) is the block Toephtz matrix, 



/ 



Bl(A) 



Ho n_i 

fli flo 

n^.i . . . 



no 



with 



1 



Ui = —J^dzz-'-'^z), 



where the matrix generator $(-2) is defined by the equations, 



and (h(z) 



iX 4>{z] 



J 



Xl {l-Xrz){l-X2Z-^) 
X,{l-Xlz-^){1-X*,z)^ 



(11) 



'121 



(13) 



z-plane with z = x + i y 




Case la&2 



Case lb 



Figure 2. The polygonal line E (oriented as it is indicated) separates the complex z 
plane into two parts: the part f2+ which lies to the left of S, and the part fl- which 
lies to the right of T,. Curve S is the unit circle in anti- clockwise direction. Cuts 
Jii J2 for the functions (j){z),w{z) are labeled by bold on the line S. The definition of 
the end points of the cuts A... depends on the case: Case la: \a = ^1, As — , 
— ^2, o-nd Xd — '^i^ ■ Case lb: Xa = Xi, Xb — A^^, Ac — ; o,i^d Xd — A2. 
Case 2; A^ = Ai, A^ = A2; Ac = Aj ^, and Xjj = Aj^ ^. ^ 

We fix the branch of 0(2;) by requiring that 0(oo) > 0. We use * to denote complex 
conjugation, and S is the unit circle shown in Fig. |21 The points Ai and A2 are defined 



differently for the different values of 7 and h. 
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There are three different cases ! 

Case la is defined by inequality 4(1 — 7^) < h? < A . 

It describes moderate magnetic field. 

Case 2 is defined by h > 2. This is strong magnetic field. 

In both cases Ai and A2 are real and given by the formulae 

h - Jh^ - 4(1 - 72) 1 + ^ 

2(1+7) ' ^-T^^" (") 

Case lb is defined by h'^ < 4(1 — 7^). 

It is weak magnetic field, including zero magnetic field. 

Both Ai and A2 are complex and given by the equations 



/i -iJ4(l -72) - 

2(1 + 7) • = 

Note that in Case 1 the poles of the function (j){z) (Eq. ITSj) coincide with the points 
Xa and A^, while in Case 2 they coincide with the points Xa and Ac. 

We are going now to formulate our main results. To this end we need to introduce 

the Jacobi theta-function, 

00 

^3(5) = J2 e"^™'+2"^"''. (16) 

n=— 00 

We remind the following characteristic properties of this function (see e.g. |39j): 

93{-s)=9s{s), 9s{s + l) = es{s) (17) 
^3(s + r)=e-'^™%(s) (18) 
+ + i + = 0, n,mGZ (19) 

The modulus parameter r, in our case, is determined by the physical quantities h and 
7 according to the following equation 

r = zro, ro = I{k')/I{k) (20) 
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where I{k) denotes the complete elhptic integral of the first kind, 



dx 







also k' = Vl — k'^-i and 



^(V2)2 + 72 - 1 / 7, Case la 



k = \ Jl - 72 - (V2)2 / Jl- (V2)2, Case lb 



7 / ^(V2)2 + 72-l, Case 2 



Define also the function 



/?(A) = — lnA±l, 
and the infinite sequence of real numbers, 

\m = tanh + ^''"0, m > 0, 

where ex = 1 in Case 1 and a = in Case 2. Observe that 



< A„ < 1, A, 



m — > 00 



:2i) 



(22) 



(23) 



in view of Eq. (fT^. the points ±Am are zeros of ^3 (/5(A) + ^ 



Theorem 1. Let f2 6e t/ie complex A - plane without arbitrary fixed neighborhoods 
of the points A = 00, A = ±1 and A = iA^, = 0,1,... . Then the Toeplitz 
determinant Dl{X) admits the following asymptotic representation, which is uniform 
in \ eVL. 



2A ^ d , 
T^^+dA^" 



^3 ( m + ^3 (/?(A) - ^ 



(24) 



-Oip 



L ^ 00. 



Here p is any positive number satisfying the inequality, 



1< p< \Xc\. 
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This theorem shows that in the large L hmit, the points ±Am (I23|) are double zeros 
of the Dl{X). More precisely, we see that in the large L limit the eigenvalues z/2m and 
z/2m+i from ©, dH) merge to A^: 

^2m, ^2m+l ^ ^m, (25) 

which in turn implies our main result: 

The limiting entropy, S{pa), of the subsystem can be identified with the infinite 
convergent series, 

OO OO r) 

S{pa)= E HiXJ= J2 (1 + A™)ln— — (26) 

m=—oo m=—oo ' 



It is worth mentioning that relation ()25|) also indicates the degeneracy of the 
spectrum of the matrix Bj, and an appearance of an extra symmetry in the large L 
limit. 

Observe that equation ()1U|) can be also rewritten as 

Sl{pa) = lini dAe(l + e,A);^in(z}^(A)(A2-l)-^) , (27) 

e^o+ Am Jr' dA ^ ^ 

because of the following identity 



2A_ 

A2 



lim i dAe(l + e,A)— ln(A2 - 1)-^ = L lim i dAe(l + e,A) — - 
€^0+ Jr' dA e^o+ Jr 1 - 

= 27riL lim ((2 + e) In + e In = 0. 

.^0+ ^ 2 2/ 

In this paper, we will show that series (|26|) coincides with the result of the following 
double limit procedure 



S{pa) = lim 

£—♦0 + 



lim^i dAe(l + e,A)Ain(D^(A)(A^-l) 
L^oo 47ri Jr' dA ^ 



-L 



X2i 



here the contour V is depicted in Fig^ This result can be alternatively written down 
as the following integral. 



2 Ji 



Entropy of XY Spin Chain and Block Toeplitz Determinants 9 

We conjecture that in fact tlie limits in (j28|) can be intercfianged. After the shorten 
version of this paper appeared in quant-ph, I.Peschel simphfied our expression for 
the entropy for non- vanishing magnetic field [Cases la and 2]. He used the approach 
of ^21- He showed that in these cases our formula is equivalent to formula (4.33) 
of [12]. Moreover, I. Peschel was able to sum it up into the following expressions for 
the entropy. 

e \ 



SipA 

in Case la, and 

S{pA 



1 

12 



In 



IQk' 



1 



k^\ 4I{k)I{k') 



TT 



+ ln 2, 



(30) 



In 



16 



TT 



(31) 



(A;2A;'2) 

in Case 2 . Here, I{k) denotes the complete elliptic integral of the first kind, 
k' = Vl — k^, and 



k 



{h/2f + 72-1/7, Case la 



(32) 



7 / ^(V2)2 + 72 - 1, Case 2 
These describes cases of moderate and strong magnetic field. 

We used the equation to calculate limiting entropy for weak magnetic field. 
In Case lb we derived: 

A;2 \ 



In 



16A;' 



1-/^2/4-72 
1 - h^/A 



\ 



iff 0</i<2Vl-7 



_ e\ Ai{k)i{k') 



+ ln 2, 



(33) 



Note that this case includes zero magnetic field. 

These expressions helped us to study the range of variation of limting entropy. 
Together with Dr. Franchini we found that the entropy has a local minimum at the 
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boundary of Cases la and lb: 

S{p) =ln2 at 1^1^+72 = 1 ^ 

At this boundary the ground state is doubly degenerated, but the rest of energy levels 
are separated by a gap. 

Note that the absolute minimum of asymptotic entropy min5'(p) = is achieved 
at infinite magnetic field corresponding to = in Case 2 [the ground state is 
ferromagnetic] 

2. The Asymptotic of Block Toeplitz Determinants, Widom's Theorems. 

Our objective is the asymptotic calculation of the block Toeplitz determinant -Dl(A) 
or, rather, its A -derivative -^InDL^X). A general asymptotic representation of the 
determinant of a block Toeplitz matrix, which generalizes the classical strong Szego 
theorem to the block matrix case, was obtained by H. Widom in [S3] (see also more 
recent work ST and references therein). Here is Widom's result. 
Let 

00 

k=—oo 

be an p X p matrix - valued function defined on the unite circle, \z\ = 1, and satisfying 
the following conditions. 

(i) 11^11 = Er=-oo ii^fcii + {Er=-oo mivk^y^' < 00 

(ii) det {p{z) 7^ 0, A||2|=i argdet = 

Consider the block Toeplitz determinant generated by (p{z), i.e. Dilf] = detTL[ip] 
where TlI^p] = {^Pj-k), j, = 0, L - 1. 
Theorem 2. {[56 )Definel 

G[^]=exp\^f lndet^(z)-|. (34) 
[ 2m J\z\=i z J 

I The integration along the unit circle is always assumed to be done in the positive, i.e. anti-clockwise 
direction. 
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Then the limit 

= lim J^lM^ , (35) 
exists. Moreover, the following general formula can be written for the quantity E[ip], 
E[ip] = det {T^[ip]T^[ip-']) . (36) 

This quite beautiful theorem is not very efficient in concrete apphcations. The 
remarkable fact though is that the proof of theorem 2 is based on an auxiliary fact 
which was established in the preceding work of Widom and which can be made an 
efficient tool for a large class of symbols ip, including the matrix function ^{z) from 
fll2|l which we are concerned with in this paper. 

Theorem 3 (theorem 4.1 of [55 ) Suppose that, in addition to the conditions of 
theorem 2, the matrix symbol ip^^ admits the Weiner-Hopf factorization, 

ip-^{z) = u+{z)u_{z) = v_{z)v+{z), (37) 

where the subscribes "+" and indicate the analyticity inside and outside of the unite 
circle, respectively. Suppose also that ^p{z) can be included into a differentiable family, 
A — s> V9(z, A). Then, \\iE[(p\ is a differentiable function of X and in fact. 



(i . , if {/ , , . , s , , , -.\ dip{z,X) 

dX 



— In E\ip] = — / trace 
dX 27r J\z\=i 



dz, (38) 



{u'^{z)u-{z) - v'_{z)v+{z)') 

where (') means the derivative with respect to z. 

Let Q be the set on the A - plane introduced in theorem 1. In the next section 
we will show that for every A G fi, the matrix-valued function $(2) admits the explicit 
Weiner-Hopf factorization, and hence the second Widom theorem - theorem 3 above, 
is applicable. Indeed, for the function ^{z) this theorem can be specified as follows. 

Theorem 4. Let X E Q then the following asymptotic representation for the 
logarithmic derivative of the determinant -Dl(A) = detTi[$] takes place. 



^ ' \K{z)U~\z) + V-\z)Vi{z))<^'\z)]dz (39) 



271 



trace 

z\=l 
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12 



+rL{\), 



where the error term ri(A) satisfies the uniform estimate, 



rLW\<Cp 



Xen, L>1, 



(40) 



and p is any positive number such that 1 < p < \Xc\- In / t,V.^) the 2x2 matrix-valued 
functions U±{z) and V±{z) solve the following Weiner-Hopf factorization problem : 



(a) U^{z) and V^{z) (Uj^^Iz) and V-^{z)) are analytic outside (inside) the unit circle 

(ill) [/_(oo) = V_(oo) = /. 

Remark . The replacement of the factorization of the inverse symbol by the 
factorization of the symbol itself has been made for a purely technical reason. 

Asymptotic representation plays an important role in our analysis as, in fact, 
its starting point. The truth of the matter is that we had derived formula (jH^ before we 
became aware of Widom's second theorem. In our original derivation of and proof 
of theorem 4 we used an alternative approach to Toeplitz determinants suggested by P. 
Deift in |31]. It is based on the Riemann-Hilbert technique of the theory of "integrable 
integral operators", which was developed in [HS], [01] for evaluation of correlation 
functions of quantum completely integrable (exactly solvable) models §. It turns out 
that, using the block matrix version of 33 suggested in jSTj, one can generalize Deift's 
scheme to the block Toeplitz matrices, which provides a rather simple derivation of 
equation ()39p together with the estimation of the error term indicated in theorem 4. 

The Riemann-Hilbert approach has already proved its usefulness in the theory of 
Toeplitz and Hankel determinants with scalar symbols (see e.g. jSZI, [Z2j, [HH], |59] . 

[UU]). Therefore, we believe that the block- version of the Riemann-Hilbert scheme is 

§ In its turn, the approach of is based on the ideas of [31] • Several principal aspects of the 
integrable operator theory, especially the ones concerning with the integrable differential systems 
appearing in random matrix theory, have been developed in )H5j . Some of the important elements of 
modern theory of integrable operators were already implicitly present in the earlier work 



(i) ^z)=U+{z)U^{z) = V^{z)V+{z) 2GS 
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worthwhile to present. Indeed, although in the relatively simple case of smooth symbols 
theorem 4, as it turned out, is a direct (up to the error term estimation) corollary of 
one of the old Widom's results, it is conceivable that in more challenging situations of 
the matrix symbols with singularities the Riemann-Hilbert scheme might become very 
useful, as it has been already the case for the scalar symbols. 

Having all the above reasons in mind, we have decided to include in the paper 
our original proof of theorem 4. We will do this in the next section providing all the 
necessary facts concerning integrable Fredholm operators and the Riemann-Hilbert 
method of their analysis. 

The explicit factorization of the matrix ^{z) will be performed in section 4 and 
the final proof of our main result - theorem 1, will be done in section 5. 

3. The Riemann-Hilbert approach to the block Toeplitz determinants. An 
alternative proof of theorem 4. 

3.1. The Fredholm Determinant Representation 

Let fj{z) and hj{z), j = 1,2, be 2 x 2 matrix functions. We introduce the class of 
integrable operators K defined on L2(H,C^) by the following equations (cf. |SZ!), 

{K X){z) = £K{z,z')X{z')dz' for X e L2, (41) 

where 



h{z) 



^ h{z) ^ 



fliz) 



(42) 



Let I2 denote the 2x2 identity matrix. Put 



h{z)=z'^l2, f2iz)=l2 (43) 

'''W--^. M.)^-^. (44) 



Entropy of XY Spin Chain and Block Toeplitz Determinants 14 

Then, essentially repeating the arguments of [HI], we have the following relation 

DLiX) = detil - K), (45) 

where / is the identity operator in L2(H,C^). The determinant in the l.h.s. is the 
2L X 2L matrix determinant, while the (Fredholm) determinant in the r.h.s. is taken 
in L2(H,C^). This relation can be proved by writing down the matrix representation 
of the integral operator / — i^' in the basis {z''ej}^oo<k<oo,j=i,2, where {61,62} is the 
canonical basic in (cf. [HI])- The matrix elements of operator I — K can be defined 
by following relations: 

{I-K)z''e^ = Y.j2%:^zUp, 0<k<L, a = {1,2^6) 

j=0 /3=1 

L-1 2 

(J - 6„ = 2' 6„ + ^ ^ <I>^7„V 6/3, k<OoTk>L, « = {1(^) 

j=0 f3=l 



Here 



With B2,(A) defined as {^p~a), < j, k < L — 1; a, f3 taking 1 or 2; {j,f3) numerating 
the rows and {k, a) numerating the columns , equation ()45|) becomes obvious. 

So we represented Dl{X) as a Fredholm determinant of the integral operator K. 
Define the resolvent operator R by 

{I-K){I + R) = I ^ R={I-K)'^K. 

Then, we have the general equation [33], [TF] . 

z — z' 

where 

= {I- K)-^f, and H = h{I - Ky^ 

where in the first equation (/ — K)~^ is understood as acting to the right, while in the 
second equation it acts to the left. 
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One of the main ingredients of the theory of integrable operators is the following 
Riemann-Hilbert representation for the functions F{z) and H{z). 

F{z) = Y^{z)f{z), zeE, (49) 

H{z) = {Y^)-\z)h{z), ZEE, (50) 

where the 4x4 matrix function Y{z) is the (unique) solution of the following Riemann- 
Hilbert problem: 

(i) Y{z) is analytic outside of the circle H. 

(ii) y(oo) = I4, where I4 denote the 4x4 identity matrix. 

(iii) Y_{z) = Y^{z)J{z), z & E where Y^{z) {Y_{z)) denote the left (right) boundary 
value of Y{z) on S (note, means: from inside of the unite circle!). The 4x4 
jump matrix J{z) is defined by the equations, 



J[z) = h + 2mf{z)h^{z) 



^ 2/2 -*^(^) -z'^ih-^'^iz))^ 



It is also worth noticing the integral formulae. 



Y(z 



Z' — Z 



and 



Y-\z)^h + l 



fiz')H^iz" 



z' — z 



-dz'. 



(51) 



(52) 



(53) 



We are going now to relate the A - derivative of det(/ — K) and the resolvent. To this 
end we first notice that 

d 



Therefore, we have that. 



dX 



K{z,z') 



1 z^(z')-^-l 



27r z — z 



^I, = -iK{z,z'){l2-Hz')) 



Entropy of XY Spin Chain and Block Toeplitz Determinants 16 

From this it follows that 



{z,z') = -iR{z,z'){l2-Hz')) 



and hence 



-^lndetTi{(/)} = -Trace {I - K)-^^K 
dX dX 



= trace (^R{z,z){l2 - ^{z))'^) dz, (54) 

where "Trace" means the trace taking in the space L2(S, C^), while "trace" is the 2x2 
matrix trace. 

Prom ()48|1 we conclude that 

, dF^(z) , 
Riz,z) = ^^Hiz), 



which in turn implies the equation 

, , ^ r ■ f fdF^(z) 

dX 



lndetTi{0} = iJjTace ('^^^^H{z){l2 - $(^))"'] dz. 



(55) 



3.3. The Asymptotic Solution of the Riemann-Hilbert Problem. 



Equation together with and (jHn|l it reduces the question to the asymptotic 
analysis of the solution Y{z) of the Riemann-Hilbert problem (1 - 3). For the latter 
we shall follow pij . 

The basic observation is that the jump matrix J{z) admits the following algebraic 
factorization, 



. where 



J{z) = M{z)Jo{z)N{z), 

V 



M{z) 



h z^[h-{^^)-\z 
Oo /: 



N{z) 



I 



\ 



I 

0. 



(l2-{^^)-\z)) h 



(56) 



(57) 



(5^ 
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and 



(59) 



Choose now a small e and define the matrix function X{z) according to the 
equations: 

X{z)=Y{z) if \z\>l + e, or |z| < 1 - e, (60) 

X{z) =Y{z)M{z) if l-e<|2|<l, (61) 

X{z) =Y{z)N'^{z) if l<|z|<l + e. (62) 

The new function has a jump accross the unit circle S with the jump matrix Jq{z) and 
two more jumps - accross the circles, 

51 : |2;| = 1 — e, jump matrix M{z) 

and 

52 : = 1 + e, jump matrix N{z). 

In other words, the original Rimeann-Hilbert problem (1 - 3) is equivalent to the 
problem 

1°. X{z) is analytic outside of the contour F = S U Si U S2. 
2°. X(oo) = I4, where J4 denote the 4x4 identity matrix. 

3°. The jumps of the function X{z) across the contour F are given by the equations 

• X_(z) =X+{z)M{z), zeEi 

• X^{z)=X+{z)N{z), zeE2 

• X.{z)=X+{z)Jo{z), zeE 

where the jump matrices M{z), N{z), and Jo{z) are defined in (jSTj), (fSHj) . and 
()59|). respectively and each circle is oriented counterclockwise. 
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Observe that the differences of the jump matrices on Si and S2 from the identity are 
exponentially small as n — > 00. This means one can expect the following asymptotic 
relation for X{z), 

X{z)r^X%z), (63) 

where the independent on L function X'^lz) solves the Riemann-Hilbert problem which 
is the same as the Y - problem but with the jump matrix Jo{z) instead of J{z). 

We notice that the function X^{z) can be found explicitly in terms of the 
2x2 matrix- valued functions U±{z) and V±{z) solving the following Weiner-Hopf 
factorization problem : 

(i) ^z) = U4z)U4z) = V4z)V+iz), ZEE 

(ii) U^{z) and VI (2;) {U+{z) and V+{z)) are analytic outside (inside) the unit circle 



(iii) f/_(oo) = ^.(oo) = /. 



Indeed we have that 

X'(Z) : 



Uliz) 

02 



02 



if \z\ < 1, 



(64) 



and 



X^iz) 



V 



if \z\ > 1. 



(65) 



In the next section, we will show that for all A G i7 the above Weiner-Hopf 
factorization can be found explicitely, in terms of the elliptic theta-functions. In this 
section though we won't need the explicit formulae for U±{z) and V±{z). The only 
information we will need in this section is the existence of U±{z) and V±{z), and the 
uniform estimate 



\UJz) 



\V+(z)\<C, Vz, VAel], 



(66) 



which is the direct consequence of the explicit formulae p27|l . p28j] . and (jHEI), (jHTjl of 
section 4. 
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We are now ready to formulate and prove the rigorous version of the formal 
asymptotic relation (jUHj) . 

Theorem 5. Let Q be the set introduced in theorem 1, i.e. the complex X-plane 
without arbitrary fixed neighborhoods of the points A = oo, A = ±1, and A = ±Am, 
m = 0,1,..., where A^ are defined in . Then, for sufficiently large L and all 

\ E VL, the Riemann-Hilbert problem (1^ - !^) has the unique solution X{z) which 
satisfies the following uniform estimate, 

x{z) {x%z)y'-h 

where p is any positive number such that 1 < p < \Xc\, and the block diagonal matrix- 
valued function X^{z) is defined by the equations (|6'^| ) and / 1 6'^) . 
Proof. Put 

R{z) ■.= X{z) {X\z)y\ 

Observe that since det Jo(-2) = 1, the function d{z) = detX°(z) solves the following 
scalar Riemann-Hilbert problem: 

• d{z) is analytic outside the unit circle S 

• d-{z) = d+{z), z 

• d{oo) = 1 

By virtue of the Liouville theorem, this implies that d{z) = detX^{z) = 1, and hence 
the matrix ratio R{z) is well defined. In fact, in terms of the matrix function R{z) the 
Riemann-Hilbert problem (1° - 3*^) can be rewritten as follows . 

1°°. R{z) is analytic outside of the contour Fq = Si U S2 (no jump across the unite 
circle S!). 

2*^°. R{oo) = I4, where I4, denote the 4x4 identity matrix. 

3°°. The jumps of the function R{z) across the contour To are given by the equations 

• R_{z) = R+{z)M{z), zeEi 

• R4z) = R+{z)N{z), zeE2 



<C / zee, Xen, L>Lo, (67) 
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where the jump matrices M{z) and N{z) are defined by the formulae 

M{z) = X\z)M{z)[X\z)]-\ N{z) = X\z)N{z)[X\z)]-\ (68) 
Let us denote Go{z) the jump matrix of the function R{z), i.e. 

[ M{z), zeEi 

Go{z) = (69) 
[ N{z), z e S2 

Then, taking into account the estimate and the expressions for matrices M{z) and 
N{z) we immediately arrive at the estimate, 

1/4-^0(^)1 < Cp"^, l<p<|Ac|, zeTo, Xen, L>1. (70) 

This in turn implies the L2 fl Loo - norm estimate: 

11/4 - G'o||L2(ro)nLoo(ro) < Cp~^, (71) 
which, by standard arguments based on the analysis of the relevant singular integral 
equation (see e.g. [SB] or Appendix D of [Z2]), yields the estimate (jHTjl . For the reader's 
convenience we will now present those arguments. 

According to the general theory of Riemann-Hilbert problems (see e.g. jHIj or [221) 
the solution R{z) of the Riemann-Hilbert problem (1*^° - 3*^°) is given by the integral 
representation, 

1 r dz' 
R{z) = /+;—/ r{z') ih - Go{z')) , z ^ To, (72) 

ZTll JTo Z — Z 

where r{z) = R-^^{z) solves the singular integral equation 

r{z)=I + C4r{h-Go)]{z), zeTo, (73) 
and C+ is the corresponding Cauchy operator: 



iG+h){z) = Urn f 



h{s) ds 



,^'e(+)-side-'ro s - z' 27rz 

We note that in the case under consideration this statement is an (almost) immediate 
consequence of the Cauchy theorem applied to the Cauchy integral. 
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in conjunction with the jump relations 3°°. 
Put 

ro{z) = r{z) - h. 

Then equation ()73|) can be rewritten as an equation in L2(ro), 

(l-K)ro = /, (74) 
where the function f{z) and the operator K are defined by the formulae, 

/ = C+(J-Go), 

K:h^C+[h{h-GQ)l heL2{To). 
The L2-boundness of the operator C+ (see e.g. and the estimate (f7T|) imply that 



U2(ro)^L2(ro) < l|C'+IU2(ro)^L2(ro)ll-^4 - GollLooCro) - ^(p 



and 

"^'|L2(ro) ^ ll'-^+llL2(ro)^L2(ro)ll-'4 - 'j-o||L2(ro) 



< l|C'+IU2(ro)--L2(ro)ll-^4 - G'o||L2(ro) - 0{p ^) 



Therefore, equation (f7I|l is uniquely solvable in L2(ro) for sufficiently large L, and its 
solution satisfies the estimate, 



ko||L2(ro) = 0{p ^) 



or, in terms of r(z) 

"'|L2(ro) 

This equation together with (f72|) allows us to estimate R{z) 



I-R(^) -h\< CWh - Col I Loo (Fo) ^. , / . -r ^ + C\\h - Col liaCFo) 1 1-^4 -r\\L2{ro)- 



'dist(2;ro) ->^M-2,.o>M-. ■ "--^^oMist (z;ro 



<C- ^ 



dist {z; Tq) 
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This leads to the announced estimate (j67|) for all z outside of some neighborhood of 
the contour Fq. Since of the flexibility in the choice of the circles Si and S2 (i.e. the 
flexibility in the choice of e in ()60I -IU ^ ). the estimate is extended to the neighborhood 
of the contour Fq. This completes the proof of the theorem. 

3.4- The Proof of Theorem 4- 

Let us observe that in terms of the function X{z) the equations (j49p and (|5(jp. which 
determine the functions F{z) and H{z), can be transformed as follows. 



F{z)=X4z)M-'iz)fiz) 



= X+{z) 
and, similarly 



( h -z^{h-{^^)-' 



V 



\ 




( 




= X4z) 


/ 


\ ^2 j 


V 



z^(^^)-Hz) ^ 



H{z) {h - Hz))-' = {Xl)-\z)M\z)h{z) {h - 
( 



This implies that 



h O2 

z^{h-^-\z)) h 



\ ( \ 

Z 12 



V 



27ri 



iXl)-\z) 



( .y-Lj \ 

Z 12 



V 



<^-Hz) 



( '^El^H{z)(h - MZ))A = —-iT<^-Hz) + —tiT^-Hz) 



tr \ —^H{z){h~Hz)) . 

\ dZ J iTil z 



2711 dz 



dz 



V 



-z^^-^{z) -^-^{z) 

and hence, our basic relation for the Toeplitz determinant ()55|) can be re-written as 



d ,^/^^ L r 1 , .dz 
— \nDLi\) = — / tT^-\z) — 
aX ZTT Js z 
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H / tr 



dXljz) 
dz 



T\-l, 



( 



V 



dz. (75) 



The second term in the r.h.s. of equation (fTSj) can be spht into the following three 
integrals, 



1 

2^ 7s 



tr 



dXl{z) 
dz 



{X. 



T\-li 



I 



\ 



-z^<l>'V 



tr 



27r 



H / tr 

271 k 



dz 



dXliz) 
dz 



V 



0, 



O2 



dz 



dz 



T\-li 



O2 O2 

-z^<^-^iz) O2 



1 



tr 



dXl{z) . I O2 z-'^h \ 



2tx 7s 

Ti + T2 + T3. 



dz 



\x\ 



0, Oo 



/ 

dz 



dz 



The integral T2 can be replaced by the integral over the contour Si 



T2 = — / tr 
27r Jsi 



and hence estimated as 



dXl{z).^^.^ 



dz 



iXi)-\z) 



O2 O2 

-z^^-'^(z) Oo 



dz. 



|T2|<Cp"-^, l<p<\Xc\, XeQ, L>Lo 



(76) 



(77) 



(7J 



(We notice that by virtue of theorem 5, the matrix function X{z) is uniformly bounded) 
Using the jump relation, X_{z) = X^{z)Jo{z), the integral T3 can be rewritten as 



dXTjz) 
dz 



I 



(X 



T\-li 



O2 z-^<!>-^{z] 

02 02 



dz. 
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and, similarly to the integral T2, further replaced by the integral over the contour S2, 



27r Jb 



tr 



dz 



>2 
O2 



0, 



dz. 



(79) 



which in turn yields the same estimate ()78|) as in the case of integral T2, 

|T3|<Cp"^, l<p<|Ac|, Xen, L>1. (80) 

Finally we notice that, by virtue of theorem 5, the matrix function X^{z) in the integral 
Ti can be replaced, within the same error (jHUj) . by the block diagonal matrix function 
X^{z) from equation ()fi4|l . Therefore, we derive from (f7H|l the following asymptotic 
representation for the logarithmic derivative of Dl{X). 



H / tr 



dz 



V 



02 



dz+rL{X){81) 



where the error term ri(A) satisfies the estimates, 

\rL{X)\ <Cp-^, l<p<|Ac|, AGfi, L>1 
Taking into account formula ()64|1 for X^{z) we can simplify (j8T| as follows. 



d 

dX 



\nDLi\) = ^ljT<^-\z) 



\( \dz 
z 



H / trace 

27T J\z\=l 



{U'4z)U-\z) + V-\z)Vi{z))<^-\z) 
+rL(A), 



dz 



(82) 



^3) 



or 



d 



2A 



+ — trace 

ZTT J\z\=l 



d\'" 1-A2 

{u'^{z)u:^\z) + v-\z)v:,{z))^-\z) 

+rL(A), 



dz 
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where we have used the relation 



tr$^^(z) 



1- A2 



In the last formulae, the symbol (') means the derivative with respect to z. This 
completes the proof of theorem 4. 



4. Wiener-Hopf Factorization of Matrix $(2; 

By explicit calculation, one can find that 

(l-A2)(T3<f-^(^)a3 = <f(^), (T3 

Hence, 



^ ^ 



-1 



^5) 



V^z) = a;Uz\z)a, (86) 
V+{z) = asU-\z)as{l-X'), A ^ ±1, (87) 

and one only needs the explicit expressions for U±{z). 

Our observation is that, for all A outside of a certain discrete subset of the interval 
[—1, 1], the solution to the auxiliary Riemann-Hilbert problem (i-iii) exists; moreover, 
the functions U±{z) can be expressed in terms of the Jacobi theta-functions. Indeed, 
the auxiliary Riemann-Hilbert problem (i-iii) can be easily reduced to a type of the 
"finite-gap" Riemann-Hilbert problems which have already appeared in the analysis of 
the integrable statistical mechanics models (see [SB])- Before we give detail expressions, 
let us first define some basic objects: 



w 



(z) = v'(---Ai)(2-A2)(2-Aj')(2-Ar'), 



where w{z) is analytic on the domain C\ { Ji U J2} shown in Fig. |21and fixed by the 
condition: w(z) — > z"^ as 
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z ^ oo. Next we define 

r = - — -, c = 2 — -, (90) 
c Jxb w[z) J\a W{Z) 

c y J\a w{zj J c Jxa w{z) 

A{z) = - / ^dz, = / du{z), (92) 

Points Ayi, As, Ac, Ad and cuts Ji, J2 and curves E and S are shown in Fig. El We shall 
also need: 



Aq = lim 



A(^)-lln(z-Ai) 



(93) 



Here, the contours of integration for c and S are taken along the left side of the cut Ji. 
The contour of integration for r is the segment [A^, Ac]. The contours of integration for 
K and in ()93p are taken along the line S to the left from A^; also in (jUH)) . arg(2 — Ai) = vr. 
The contours of integration in the integrals A{z) and uj{z) are taking according to the 
rule: The contour lies entirely in the domain fi+ (^-) for z belonging to (1^-). It 
also worth noticing that, for the general reasons (see e.g. [ZOj; see also Chapter 2 of 
[7T]). we always have that Rezr < 0. Moreover, as we will see later on, the following 
inequality takes place 

ir < 0. (94) 

It also should be noted that uj{z) is an elliptic integral of the first kind, normalized by 
the period condition 

2 [^^ duj{z) = 1. (95) 
Jxa 



and A{z) is an elliptic integral of the third kind, normalized by the conditions 

IXa 



2 T'' dA{z) = -TTi, (96) 

JXa 



and 



A(z) ~ -Inz for z ^ 00. (97) 
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The rules formulated above define the integrals A{z) and uj{z) as the single- valued 
analytic functions in the domains f2+ and Let us denote and S/j the parts of the 
contour S which lie to the left of = Ai and to the right of \d = ^i^, respectively || . 
Then, on the polygonal line S the integrals A{z) and uj{z) have jumps which can be 
described by the equations, 

uj+{z) - uj^{z) = and A+(^) - A_{z) =0 for z G S^, (98) 

uj+{z) + uj_{z) = and A+(z) + A_{z) = for ;z G (A^, A^) (99) 

uj+{z) - uj^{z) = 1 and A+{z) - A_{z) = -ni foi ze{\B,Xc) (100) 

uj+{z) + uj^{z) = T and A+{z) + A^{z) = 2 / dA{z) = -2^ / duj{z) 

for z G (Ac, Ad) and integration is taken over the ray (101) 

uj+{z) - uj-{z) = and A+{z) - A^{z) = -Txi for z eT^o (102) 

The only relation above which is not straightforward is the last equation in (jlOll) . It 
easy follows, however, from the application of the classical Riemann bilinear relations 
(see e.g. [ZOj o^' [H]) to the periods abelian (elliptic) integrals uj{z) and A{z). 
Let 



oo 



e^{s) = e"'™ (103) 

n=— oo 

be the third Jacobi theta-function. We remind again the basic properties of the 
function6'3(s) (see e.g. P^): 

e,{-s) = e^{s), e;{s + i) = e^{s) (io4) 

e^{s + t)= e-^'^-^^''e-i{s) (105) 

fn + mr + ^ + ^) = 0, n.me'L (106) 



2 2 

We also introduce the 2x2 matrix valued function Q{z) with the entries 



en(^) = (^-Ai)-^e^(^) 



II In particular, for the Cases la and 2, we have that and Ejj are the half lines {^oo,Xa) and 
(Ad,+oo), respectively. 
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X 



03 

-(z - Ai 

03 {00{Z) 


[u{z) + f ) 

-/9(A) + K 


err \ 
2 y 


03 

-(z - 
6*^ (uj(z) 


[uiz) - 

1 \ J ' 


err \ 
2 


^3 

(^-Ai)- 

^3 {UJ{Z) 


[uiz) - f ) 
-/3(A) 


+ f) 


03 


(^(^) + f ) 





(107) 

where a = 1 in Case 1 and a = in Case 2, and /5(A), uj{z) and /t are defined 
in Eqs. fl88H93p . The branch of (z — Xi)~^ is defined on the 2;-plane cut along the 
part of the hne S which is to the right of Ai = Aa, and it is fixed by the condition 
arg(2; — Ai) = tt, if 2; — Ai < 0. 

The matrix function Q{z) is defined on C\ S. However, taking into account 
equations (jHEl - ll()2j) describing the jumps of the integrals lj{z) and A (2) over the 
line S and the properties p04p and p05|) of the theta function, one can see that 6(z) 
is actually extended to the analytic function defined on C\ { Ji U J2} which satisfies the 
jump relations 

e+(z) = e(z)ai 



A = i 



A+ 1 
A - 1 



z e Ji 




(108) 


' zeJ2. 




(109) 


\ 


(0 1] 




, 0-1 = 




(110) 


/ 







The end points of the intervals Ji,2 are singular points of the function Q{z). Specifically, 
one can notice that 



e(z) = O 



.Vz-Xi^ 



as z ~ Ai, 
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as -2 ~ A2 ^■ 



The last estimate follows from the fact that, according to ()95p and ()106|) . the only 
zero of the denominators of the right hand sides in ()l()7j) is the point 2; = A2 ^ (which 
is Xb in the Case 1, i.e. when a = 1, and Ac in the Case 2, i.e. when a = 0). At 
the end points z = X^^ and z = X2 the function Q{z) is holomrphic with respect to 



the variables y z — X-^^ and y/z — X2, respectively, although det &{z) vanishes at these 
points. In more detail, the following representations take place at the points z = Xa, 
Xb, Ac, and Ad : 



in a neighborhood of z = Xa, cut along Ji 



Qiz) = eA{z)iz - Xa)^ 



1/2 




1 1 



in a neighborhood of 2; = A^, cut along Ji, 

/ 



Q{z) = eB{z){z-XB)\ 



1/2 -a 




(i -i\ 
1 1 



(112) 



in a neighborhood of z = Ac, cut along J2 

( 



Q{z) = Qc{z){z~Xc)\ 



a -1/2 




/A 



J 

(1 -A 
1 1 



A-i 



(113) 



in a neighborhood of 2; = A^, cut along J2 

/ 



Q{z) = Qo{z)iz-XD)^ 



1/2 




J A 



1 -l\ 



1 1 



A-\ 



(114) 



In these formulae, the left matrix multipliers Qa{z), Qb{z), Qc{z), and Qd{z) are 
matrix - valued functions holomorphic (with respect to the variable z itself) and 
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invertible as matrices at the points z = \a, z = Xb, z = Xc, and z = Xd, respectively. 
To prove, say, representation let us introduce the model function. 



^ -1/2 ^ 



eo(^) := {z - Xa) V 







1 1 



This function is well defined in a small disk centered at and cut along Ji. By 
a simple direct computation we observe that the function Qo{z) satisfies accross the 
portion of Ji lying in the disk exactly the same jump condition as the one indicated in 
(ll(J8j) . Therefore, the matrix ratio, 

QAiz) ■.= e{z)e^\z), 



does not have any jumps at all in the neighborhood of A^. Hence, z = A^ is its isolated 
singular point. From the explicit formulae ()107|) it follows that, at worst. 



\Qa(z)\ < 



C 



\z-Xa\'/''' 



which implies that z = Xa is in fact a removable singularity. Equation pilj) follows . 
The other equations from the list (jllll - ITTH) can be proven in the similar way. 

Strictly speaking, we have not completed the proof of equations pill - (TT^ . We 
need to show that all the left multiplies are invertible matrices. To see this we note 
that 

O3 



enfoo) 



+ f) 



©22(00) 
©12(00) 



03 + f ) 
,A„ ^3(/3(A)-f 

^^3 («: + f ) 
©21(00) = , 



(115) 

(116) 
(117) 
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To prove the latter equation, consider the ratio 
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:ii8i 



ip{z) := det e{z)/(l){z). 

The jump conditions ()in8|l and (jlOQj) imply that the function iIj{z) has no jumps while 
the representations ITT^ show that neither it has poles. Hence, this function is 

a constant and we have that 

ip{z) = tpi^oo) = det 9(oo)w — . 

V Ai 

This proves pi8|) . In its turn, equation ()118p implies the invertabihty of Qki^k) for all 
k = A, B,C, D. In addition we also have that 

det Q{z) ^ 0, for all z i {X[^ , A2}. 



In addition to matrix- valued function Q{z), we introduce the matrix- valued 
function 



Q{z) 



V 



>[Z) -(p[z) 
i i 

Note that Q{z) diagonalizes original jump matrix $(z): 

^z) = Q{z)AQ-\z) 
and Q{z) is analytic on C\ { Ji U J2} and 

Q+{z) = Q_{z)cri, zeJiUJ2. 



(119) 



(120) 



Moreover, at the end points of the intervals Ji^2, which are the only singularities of 
the matrix function Q{z), it behaves in the way similar to the behavior of the function 
e{z), i.e. we have (cf. (ITTT]- ITTi|l ). 
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Q{z) = Qa{z){z-Xa)^ 



1/2 




1 1 



in a neighborhood of z = \b, cut along Ji, 

/ 



Q{z) = Qb{z){z-Xb)\ 



1/2 -a 




(122) 



^ -1^ 
1 1 



(123) 



in a neighborhood of z = Ac, cut along J2, 

I 



a -1/2 




1 1 



(124) 



in a neighborhood of 2; = Ad, cut along J2, 



^1/2 ^ 



Q{z) = Qn{z){z-\n)^ ° ^ I 



(125) 



with some holomorphic and invertible at the respective points matrices Qa{z), Qb{z), 
Qc{z), and Qd{z). We are now ready to present the solution U±{z) of the Riemann- 
Hilbert problem (i-iii). Put 

A = Q{oo)A-^e-\oo), (126) 



and define 



U4z) := Ae{z)AQ-\z), \z\ > 1 
U+{z) := Q{z)&-\z)A-\ \z\ < 1. 



(127) 
(128) 



By virtue of Eq. (|12(jp . we only need to be sure that U^{z) and U^{z) are analytic and 
matrix invertible for \z\ > 1 and \z\ < 1 respectively. From the jump properties p08l 
- ll(J9j) of 6(z) and (jl21|) of Q{z) it follows that U± have no jumps across Ji^2, and 
hence they might have only possible isolated singularities at Ai_2, ^il- The analyticity 
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at these points follows immediately from the representations (jllll - I114j) and (jl22l - 
I125p which show that the singularities, which the functions Q{z) and Q{z) do have at 
the end points of the segments Ji^2, are canceled out in the products ()127|) - p28|) . 

The excluded values of A for which the above construction fails are A = ±1 and, 



in view of Eq. (jllSp , the zeros of ^3 (/5(A) + i.e. (see (HUSD) 
± Am, Am = tanh (m H — ^ tttq, m > 0, 



(129) 



where, 



To = -tT 



f Ac dz 

.J\b w(z) 
I ^-^ 

r^B dz 

J^A w{z) 



> 0. 



We conclude this section by noticing that, using the standard reduction (see e.g. 
j) of an arbitrary elliptic integral to the canonical elliptic integrals, we can rewrite 
the above expression for tq as 

To = I{k')/I{k). (130) 

Here I{k) denotes the complete elliptic integral of the first kind. 



m 



dx 



^{l-x^){l-k^x^) 



k' = Vl — k^, and 



(/i/2)2 + 72 _ 1 / Case la 



k = < Jl - 72 - (/i/2)2 / Jl - {h/2)\ Case lb 



7/^(V2)2 + 7'-l, Case 2 



5. The Proof of Theorem 1. Evaluation of the Entropy. 



:i3ii 



Denote 



six) := — / trace 
^ ' 27ri|.|=i 



{y'^{z)U-\z) + V~\z)Vi{z)) ^^\z)\ dz (132) 



The proof of the theorem 1 will be achieved by showing that 

d 



s{\) 



dX 



In 



(133) 
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From equation (|H7j) we obtain 



Denote 



s(A) = — / trace \U'^{z)U-\z) (<l>'\z) - cx3<^-\z)a3)] dz. (134) 
Itx J\z\=i '- ^ ^ 



m{z) := ^-\z) - a,^-\z)a3 



-<j){z] 



V r\z) j 



From equation (|128p . we have 



U~\z) = AQ{z)Q-\z), U'Az) = Q\z)Q-\z)A-' + Q{z){Q-^)\z)A-^ 



and from equation ()119|) . we have 



( A.-l( 



Then 



s(A) = — / trace U'{z)UT^{z)miz] 
27r J\z\=i L 

d 



dz 



trace 



z\=l 



Q-\z)-Q{z)a, 
dz 



dz. 



Let 



aiz) := trace 



d 



e-(.)-e(.).3 



(135) 



(136) 



Then, one can easily see that the jumps of a{z) across Ji U J2 satisfy the equation, 

Oi+{z) = —a^{z). (137) 

From the exphcit formulae ()107|) for the B-matrix we derive the following asymptotic 
representation (actually, the Laurent series) at z = 00 



Biz) = e(oo) / 



(138) 



where G(oo) is the diagonal matrix defined in pi5|) - ()115p and the matrix coefficient 
is given by the equations. 



Mil = Ai + + ^ ^ 



2 



(139) 
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M22 = Ai + — + ) f + y— 140 

2 ce,[^ + ^) ce,[/3{x)-^-^ 



I,, = _e-2Ao^A^_^ ^. (142) 



12 — ~e 7 r 7 r- IJ-^J-J 

^3(/3(A) + f) ^3(«:-f) 
^^ ^3 (/3(A) + 2K-f) ^3 + 

In the above equations, Ai is the coefficient of term of the expansion of the integral 
A (2;) at z = 00 and 6*3 means the derivative of the theta function with respect to its 
own argument (not the derivative with respect to A !), i.e. 

In its turn, representation ()138p immediately implies the following asymptotic behavior 

of a{z) as z — > 00, 

aW = ^ + 0(l), .-00. (143) 

where 

m(A) = —trace [Mcr3] 

(TT 

2 



c 



.^3(/3(A) + f) ^3 (/9(A)- f) 



The representations (jllll) - pi4|) of the function G(z) at the end points show that at 
the end points, i.e. at the points Ai, A2, A|f \ A2 ^, the function a{z) has at most square 
root singularities. This, together with ()137p and ()143p yield the final explicit formula 
for a{z) (cf. the derivation of equation ()118|) given above): 

a(.) = (145) 

W[Z) 

Equation ()145|1 means that equation ()135j) can be rewritten as 

2i m(X) r^B dz , , , 

^ A = / 146 

TT 1 — A^ Jxa w[z) 
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where the integration is taken on the left side of the cut [A^, A^]- By virtue of (jl44p 
and formula for the constant c, the last equation can be transformed into 



.(A) 



TtI- A2 



err 

2 



»3 /3(A) + f »3 ,fl(A) " 



Taking into account that 



we arrive to the formula ()133|) . Theorem 1 is proven. 



We are now going to show that the two expressions for the entropy S{pa), i.e. 
equations and (PH|) coincide. We will proceed by defining S{pa) by equation (PS|) 
and showing that formula (pUj) holds. To this end we rewrite as 

S{pa) = limS^ipA), (148) 

S^ipA) = lim ^ / dA e(l + e, A);^ In (z^l(A)(A2 - 1)"^) , 
L^oo 47ri Jr' dA ^ 

and observe that, by virtue of theorem 1, 
^ £ dA e(l + e, A) A In (d,(A)(A^ - 1)"^) = ± £ dA e(l + e, A).(A) + 0{p-% 

where s(A) is given by ()133p . Therefore, 

SeiPA) = / dA e(l + e, A)s(A). (149) 
47n Jt' 

The function s(A) is analytic on the A - plane minus the points ±1 and ±Am,- Moreover, 
it is an odd function with the zero at A = oo. Therefore, its Laurent series at infinity 
is of the form, 

C3 C4 , 
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This implies that (jl49|) can be transformed as follows ( see Fig. Q). 



SeiPA) = ^<f dAe(l + e,A)s(A) 
47ri JT' 

= lim / dAe(l + e,A)s(A) 
R^oc Amir 



1 /--i-^ 
Am J-oc 



1 + e + A 



In 



1+e + A 


. \ 1 + e + A / 
+ .vrj+ 2 [in 


1 + e + A 


— iiT^ 


2 




2 





s{\)d\ 



Ani J I 



1 + e- A 



In 



1 + e- A 



1 + e- A 



In 



1 + e- A 



+ in 



s{\)dX 



= -\ I ^ \l + e + \)s{\)dX + \ r {1 + e- X)s{X)dX 

4 J-oo 4 Jl+e 

= - /°°(l + e-A)s(A)c/A, (150) 

2 Jl+e 

where in the last equation we have again used the oddness of the function s{X). 
Representing the function s(A) as 

, ^^3 (/?(A) + f ) ^3 (/?(A) - f ) 
■5(A) = ^In — ^ — -, 

dx ei (f ) 

we can perform in ()15()|1 integration by parts. Indeed, the function 

t(A) = In- 



(151) 



^3 


[m + f ) 


^3 


[m - f ) 


01 







is smooth for all A G [1 + e, +oo), and 



t(A) = 0(-), A^oo. 



This yields the following, alternative to 
S,{pa)- 

1 O3 
SeiPA) = 7: In — 

2 Jl+e 



integral representation to the entropy 



+ f)Os 



(TT 

2 



dA, 



(152) 
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It is an easy corollary of the second periodicity property (jl8|) of the theta function 

that 



Observe that 



and hence 



Also, 



In 93{ia + b) = —a^ (^1 + O ) ) , a ±00, 

dist + 6; |i + ^ + n + mr }> ) > (5o > 0. 



/?(A) G zE (modZ) , VA G (-00, -1) U (1, 00) 



dist (^/?(A)±^;|l + ^ + n + mr}) > ^ > 0. 



i/?(A) ^ ±cx), 



(153) 



as A ^ ±1 ±0. Therefore, (|153|) is applicable and we see that the function t(A) satisfies 
the estimate, 



t{X) 



2'ktq 



A-1 



A + 1 



, A^±1±0, (154) 

and hence is integrable at A = ±1. This means, we can take limit e — ^> in p52|) at 
arrive to the final integral formula for the entropy S{pa) (cf. equation (|^^) 



In 



^3 


[m + f ) 


^3 


^/^(^) - f ) 


01 (f ) 





d\. 



(155) 



To obtain the infinite sum representation (j26|) we introduce the positively oriented 
contour Tm,R in Fig El 

Then, by residue theorem. 



m 1 „ 

e{l,\k) = — f dAe(l,A)s(A) 
.1 _^ 47ri 7r„ 



fe=— m— (1— (t) 



lim — <i dAe(l,A)s(A). 
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Figure 3. Contours T^.r with small circles ±Cm around the points ±1 with the 
radius r„i = 1 — ^j^-^^X^^ + j^:^¥T^Xm+ij o,Tid the big circle C'ji centered at zero 
with radius R, the two copies of the interval [1 + rm-,R] - each lying on one of the 
sides of the cut from 1 to +oo, and two similar copies of the symmetric interval, 



The last integral can be transformed (cf. derivation of (jlSUp ) to the expression 



— / dAe(l,A)s(A) + — / dAe(l,A)s(A) 
4:711 J Cm 47ri J -Cm 



1 /■- 1-*-" 



\ roc 



1 + A 



In 



1 + A 



+ + — — In 



1 + A 



ITT 



4:TCi Jl + 





1 - A 




1-A 












+ «7r^ 




2 




2 





s{X)dX 
s{X)dX 



— (f dAe(l,A)s(A) + — / dAe(l,A)s(A) 
47ri Jcm 47ri J~Cm 

.-/ (i + X)s{X)dX + - il-X)s{X)dX 



dAe(l,A)s(A)+-^/ dAe(l,A)s(A)+^ r (l-A)s(A)t;A, (156) 

Am JCm Ani J-Cm 2 Jl + rm 
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By a straightforward calculation we find that uniformly in A G Cm the following 
estimates hold, 

/3(A) = -|--r(^m+l-0+O(e-2-o), A = 1 + r^e^'', m ^ oo, 
and hence 

dist (/5(A) ±^;|l + I+ n + mV}) > ^ > 0, 

for all A G ±Cm and sufficiently large m. This again implies the applicability of 
and, as a consequence, the estimate 

|s(A)|<— In— , Ag±C„. 

The last estimate in its turn means that the integrals over ±Cm in p56|) vanish as 
m — * oo while the integral over [1 H-r^, oo] becomes the integral over [1, oo]. Therefore 
we arrive at the relation, 

oo oo oo r) 

e(l,A,)^ ^ H{\m)^ E (1 + AJln— -— 

A;=— oo r?i=— oo m=— oo ™ 

2 ii ^ ^ ^ 2 A 

which completes our evaluation of the entropy. 

6. Some critical cases. 



The entropy has singularities at phase transitions. When r — we can use Landen 
transform (see [HH]) to get the following estimate of the theta-function for small r and 
pure imaginary s: 
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Now the leading term in the expression for the entropy (|155p can be replaced by 

S{pA)='£ + 0i^-^j forr^O. (157) 

Let us consider two physical situations corresponding to small r depending on the case 
defined on the page 2: 

(i) Critical magnetic field: 7 7^ and h ^ 2. 

This is included in our Case la and Case 2, when h > 2y/l — 7^. As h ^ 2 the 
end points of the cuts Ab — > Ac, so r given by Eq. (IHUjl simplifies and we obtain 
from Eq. ()157|) that the entropy is very large: 

^(p^) = -iln|2-/i| + -ln47^+cx) (158) 
6 3 

h — y2 and 7^0. 

Next correction is 0{\2 — /i| In^ |2 — h\). This limit agrees with predictions of 
conformal approach The first term in the right hand side of p58|) can be 

represented as (1/6) In^, this confirms a conjecture of ^21- The correlation length 
^ was evaluated in Hlj . 

(ii) An approach to XX model: 7—^0 and h < 2: It is included in Case lb, when 
Q < h < 2\/\ — 7^. Now Ab — > Ac and A^ Ad, we can calculate r explicitly. 
The entropy increases without a bound: 

S\pa) = -^ln7 + ^ln(4 -h^)^\\^2^ +00, 

7 — ^ at h<2 (159) 

correction is 0(7ln^7). This agrees with [10] • 

It is interesting to compare this critical behavior to Lipkin-Meshkov-Glick model. 
It is similar to XY model but each pair of spins interact with equal force, one can 
say that it is a model on a complete graph. The critical behavior in Lipkin-Meshkov- 
Glick was described in it is similar to XY ^ but actual critical exponents are 
different. Lipkin-Meshkov-Glick model displays interesting behavior of the entropy. 
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Both the concurrence and the entropy were studied numerically and analiticaly [in the 
thermodynamical limit] in j54j . 

Remark . The zeros satisfy an estimate: 

\Xm+i — Ami < 47rro with tq = —it. 

This means that (Am+i — Am) ^ as r ^ for every m. This is useful for 
understanding of large L limit of the XX case corresponding to 7 0, as considered 
in \40l - The estimate explains why in the XX case the singularities of the logarithmic 
derivative of the Toeplitz determinant d\nDL{X)/dX form a cut along the interval 
[—1, 1], while in the XY case it has a discrete set of poles at points ±Am of Eq. M2y\jl . 

It would be interesting to generalize our approach to the a new class of quantum 
spin chains introduced recently by J. Keating and F. Mezzadri, while study matrix 
models [Hj. 
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